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ABSTRACT 

Using a catalog of 147,986 galaxy redshifts and fluxes from the Sloan Digital Sky Survey 
(SDSS) we measure the galaxy luminosity density at z = 0.1 in five optical bandpasses cor- 
responding to the SDSS bandpasses shifted to match their restframe shape at z = 0.1. We 
denote the bands 01 u, 01 g, ° V, °\ 01 z, with A off = [3216, 4240, 5595, 6792, 8111 A] respec- 
tively. To estimate the luminosity function, we use a maximum likelihood method which allows 
for a general form for the shape of the luminosity function, simple luminosity and number evo- 
lution, incorporates the flux uncertainties, and accounts for the flux limits of the survey. We 
find luminosity densities at z = 0.1 expressed in absolute AB magnitudes in a Mpc 3 to be 
[-14. 10 ±0.15, -15. 18 ±0.03, -15.90 ±0.03, -16. 24 ±0.03, -16. 56 ±0.02] in [ 01 u, 01 g, ° V, °\ 
°' 1 z], respectively, for a cosmological model with CIq = 0.3, Oa = 0.7, and h = 1, and using 
SDSS Petrosian magnitudes. Similar results arc obtained using Scrsic model magnitudes, sug- 
gesting that flux from outside the Petrosian apertures is not a major correction. In the V 
band, the best fit Schechter function to our results has <\> % = (1.49 ± 0.04) x 10~ 2 h 3 Mpc -3 , 
M* — 5 log 10 h = —20.44 ± 0.01, and a = —1.05 ± 0.01. In solar luminosities, the luminosity den- 
sity in ° V is (1.84 ±0.04) h 10 s Lo.^ q Mpc -3 . Our results are consistent with other estimates 
of the luminosity density, from the Two-degree Field Galaxy Redshift Survey and the Millenium 
Galaxy Catalog. They represent a substantial change (~ 0.5 mag) from earlier SDSS luminosity 
density results based on commissioning data, almost entirely because of the inclusion of evolution 
in the luminosity function model. 

Subject headings: galaxies: statistics 

1. Motivation 

Establishing the low-redshift galaxy luminosity density of the universe is a fundamental measurement 
of the contents of the local universe. The last two decades, beginning with the Center for Astrophysics 
redshift survey (Huchra et al. 1983), have seen steady progress in understanding the total galaxy luminosity 
density emitted in the universe. Part of the progress has been due to measuring larger and larger numbers 
of galaxy fluxes and redshifts, from only a few hundred redshifts before 1980 to a few hundred thousand 
at present. However, equally importantly, the determination of galaxy fluxes has been steadily improving. 
The luminosity functions calculated from the CfA survey were based on "Zwicky" magnitudes: essentially, 
they were determined by visual inspection of photographic material. Automated processing of photographic 
plates provided an improved way of measuring flux; the Automatic Plate Measurement survey (Maddox et al. 
1990), from which targets were selected for the Two-degree Field Galaxy Redshift Survey (2dFGRS; Collcss 
et al. 2001), is the latest example of this method. However, deep CCD imaging, such as that performed by 
the Sloan Digital Sky Survey (SDSS; York et al. 2000), which provides a higher dynamic range and more 
linear response than do photographic plates, is yielding the highest signal-to-noise ratio measures of flux and 
the greatest surface brightness sensitivity in a large survey to date. 

A preliminary estimate of the SDSS galaxy luminosity function was performed by Blanton et al. (2001), 
using a small sample of commissioning data. At the time of writing, the SDSS has obtained photometry and 
redshifts for more than ten times the number of galaxies in the commissioning data used by Blanton et al. 
(2001); furthermore, the photometric calibration procedures have improved since that time. In addition, wc 
have developed a new and more self-consistent method for if-correcting galaxies to a fixed frame bandpass 
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(Blanton et al. 2002a). Finally, we have developed a maximum likelihood technique for fitting the galaxy 
luminosity function which allows for generic luminosity function shapes as well as luminosity and number 
density evolution (Blanton et al. 2002c). Thus, it is now time to reevaluate the galaxy luminosity function 
in the SDSS. 

Throughout this paper we assume a Friedmann- Robertson- Walker cosmological world model with matter 
density fl = 0.3, vacuum pressure Q\ = 0.7, and Hubble constant H Q = 100 h km s _1 Mpc -1 with h = 1, 
unless otherwise specified. 

In Section 2, we describe the SDSS galaxy catalog data. In Section 3, we briefly describe our method of 
fitting the luminosity function (a fuller description of the method can be found in Blanton et al. 2002c). In 
Section 4, we present our results for the galaxy luminosity density at z = 0.1. In Section 5, we test whether 
our results are sensitive to how we determine the galaxy magnitudes. In Section 6, we compare our results 
to previous work. We conclude and discuss plans for future work in Section 7. 

2. SDSS Spectroscopic Galaxy Catalog 

The SDSS (York et al. 2000) is producing imaging and spectroscopic surveys over n steradians in the 
Northern Galactic Cap. A dedicated 2.5m telescope (Siegmund et al., in preparation) at Apache Point 
Observatory, Sunspot, New Mexico, images the sky in five bands between 3000 and 10000 A (u, g, r, i, z; 
Fukugita et al. 1996) using a drift-scanning, mosaic CCD camera (Gunn et al. 1998), detecting objects to 
a flux limit of r ~ 22.5 mags. A major goal of the survey is to spectroscopically observe 900,000 galaxies 
(down to r lim w 17.77 mags), 100,000 Luminous Red Galaxies (Eisenstcin et al. 2001), and 100,000 QSOs 
(Fan 1999, Richards et al. 2002) selected from the imaging data. This spectroscopic follow up uses two 
digital spectrographs on the same telescope as the imaging camera. Many of the details of the galaxy survey 
are described in the galaxy target selection paper (Strauss et al. 2002). Other aspects of the survey are 
described in the Early Data Release paper (EDR; Stoughton et al. 2002). The survey has begun in earnest, 
and has so far obtained about 30% of its intended data. 



2.1. SDSS Observational Analysis 

The SDSS images are reduced and catalogs are produced by the SDSS pipeline photo, which measures 
the sky background and the local seeing conditions, and detects and measures objects. The astrometric 
calibration is performed by an automatic pipeline which obtains absolute positions to better than 0.1 arcsec 
(Pier et al. 2002). The magnitudes are calibrated to a standard star network (Smith et al. 2002) approx- 
imately in the AB system. There are small differences between the system output by the SDSS pipelines 
and a true AB system, amounting to Am = —0.042, 0.036, 0.015, 0.013, —0.002 in the u, g, r, i, and z bands. 
Because these effects were discovered at a relatively late date in the preparation of this manuscript, we have 
not self-consistently included these shifts in our results. Instead we have applied them a posteriori to our 
results in the °-\ 01 g, ° V, °-H, and 0A z bands. 

Object fluxes are determined several different ways by photo, as described in Stoughton et al. (2002). 
The primary measure of flux used for galaxies is the SDSS Petrosian magnitude, a modified version of the 
quantity proposed by Petrosian (1976). The essential feature of Petrosian magnitudes is that in the absence 
of seeing they measure a constant fraction of a given galaxy's light regardless of distance (or size). More 
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specifically, we define the "Petrosian ratio" IZp at an angular radius 9 from the center of an object to be the 
ratio of the local surface brightness averaged over an annulus at r to the mean surface brightness within 9: 



where 1(9) is the azimuthally averaged surface brightness profile and a\ < 1, a^i > 1 define the annulus. 
The SDSS has adopted a\ Q = 0.8 and othi = 1-25. The Petrosian radius 9p is the radius at which the TZp 
falls below a thrcshhold value TZp,u m , set to 0.2 for the SDSS. The Petrosian flux is defined as the flux within 
a circular aperture with a radius equal to Np9p, where Np = 2 for the SDSS. Petrosian magnitudes are 
described in greater detail by Blanton et al. (2001) and Strauss et al. (2002). 

Another important measure of flux for galaxies is the SDSS model magnitude, which is an estimate of 
the magnitude using the better of a de Vaucouleurs and an exponential fit to the image, photo also measures 
the flux in each object using the local PSF is a template, which is the highest signal-to-noise measurement 
of flux for point sources. Finally, photo outputs an azimuthally averaged radial profile for each object. 

For the purposes of this paper, we have implemented one more measure of galaxy magnitude: a Scrsic 
magnitude, following Sersic (1968). We fit an axisymmetric Sersic profile to the azimuthally-averaged radial 
profile of the galaxy: 



where A, r , and n are free parameters quantifying the amplitude, size, and shape of the surface brightness 
profile quantitatively. We do so by convolving I(r) with the double Gaussian fit to the seeing output 
by photo, and minimizing \ 2 with respect to the observed radial profile (using the photo catalog output 
describing the radial profile and its uncertainties, as described in Strauss et al. 2002). In practice, we fit for 
r and n in the i-band and only fit for A in the other bands, fixing r and n to their values in the i-band 
(in analogy to the model magnitudes output by photo). As found by Blanton et al. (2002b), galaxies have 
best fit Sersic indices ranging from around exponential profiles of n = 1 (which galaxies tend to be blue, low 
luminosity, and low surface brightness) to quite concentrated galaxies with n = 4-5 (which galaxies tend to 
be red, high luminosity, and high surface brightness). A de Vaucouleurs profile, which is generally ascribed 
to elliptical galaxies, has n = 4. 

While SDSS Petrosian magnitudes contain over 99% of the flux within an exponential profile, they 
contain only around 80% of the flux within a de Vaucouleurs profile (in the absence of seeing). We can 
evaluate the total flux contained in our Scrsic model; this estimate is corrected for seeing and is an attempt 
to extrapolate the profile to infinity. Comparing the SDSS Petrosian magnitudes to the Sersic magnitudes 
yields a guess of how much luminosity density we are missing because of the finite Petrosian aperture size. 
We have chosen not to use the SDSS model magnitudes for this purpose because for the versions of photo 
used for this set of photometry, the outer radius which photo considered in its model fits was too small to 
accurately model large objects. 

To supply targets selected by the imaging program for the concurrent spectroscopic program, periodically 
a "target chunk" of imaging data is processed, calibrated, and has targets selected. These target chunks 
never overlap, so that once a set of targets are defined in a particular region of sky, it never changes in 
that region. Thus, the task of determining the selection limits used in any region reduces to tracking how 
the target chunks cover the sky, which is done by an internal SDSS operational database. Within each 
target chunk, a target selection pipeline determines which objects are QSO targets, galaxy targets, and/or 
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star targets, depending on the properties of each object. The pipeline selects the Main Sample galaxies 
used in this paper according to an algorithm described and tested in Strauss et al. (2002), which has three 
most important steps: star/galaxy separation, the flux limit, and the surface brightness limit. Expressed 
quantitatively, these criteria are 

?"PSF — r model > s limit 

^pctro < nimit, and 

M50 < ^50, limit, (3) 

where r pctr o is the dereddened Petrosian magnitude in the r band (using the dust maps of Schlegel et al. 
1998), r mo dei is the model magnitude, rpsF is the PSF magnitude, and /X50 is the Petrosian half-light surface 
brightness of the object in the r-band. In practice, the values of the target selection parameters vary across 
the survey in a well-understood way, but for the bulk of the area, they are: si; m it = 0.3, riimit = 17.77, and 
A*50, limit = 24.5. As described in Strauss et al. (2002), there are many more details in galaxy target selection 
which we do not have space to discuss here. We note here that objects near the spectroscopic flux limit are 
nearly five magnitudes brighter than the photometric limit; that is, the fluxes are measured at signal-to-noisc 
of a few hundred. 

To drill spectroscopic plates which have fibers on these targets, we define "tiling chunks" which in 
principle can contain numerous target chunks (and parts of target chunks). An automatic tiling pipeline 
(Blanton et al. 2002b) is then run in order to position tiles and assign fibers to them, after which plates are 
designed (that is, extra fibers are assigned to possibly interesting targets and calibration fibers are designed) 
and then drilled. In general, these tiling chunks will overlap because we want the chance to assign fibers to 
targets which may have been in adjacent, earlier tiling chunks but were not assigned a fiber. For a target to 
be covered by a particular tile, it must be in the same tiling chunk as that tile and be within 1.49° of the tile 
center (because the edges of tiles can extend beyond the tiling chunk boundaries, a particular direction can 
be within 1.49° of the tile center but not "covered" by it as defined here). To calculate the survey window 
function, we then divide the survey into a number of small regions known as "sectors" which are regions 
which are covered by a unique set of tiles (the same as the "overlap regions" defined in Blanton et al. 2001). 
These sectors are described in samplelO; the sampling rates are calculated on a sector-by-sector basis. 

The targets are observed using a 640 fiber spectrograph on the same telescope as the imaging camera. 
The results of the spectroscopic observations are treated as follows. We extract one-dimensional spectra from 
the two-dimensional images using a pipeline (specBS v4_8) created specifically for the SDSS instrumentation 
(Schlegel et al. in preparation), which also fits for the redshift of each spectrum. The official SDSS redshifts 
are obtained from a different pipeline (SubbaRao et al. in preparation). The two independent versions 
provide a consistency check on the redshift determination. The results of the two pipelines agree for over 
99% of the Main Sample galaxies. 

As of April 2002, the SDSS had imaged and targeted 2,873 deg 2 of sky and taken spectra of approximately 
350,000 objects over ~ 2, 000 deg 2 of that area. From these results, we created a well-defined sample for 
calculating large-scale structure and galaxy property statistics, known as Large-Scale Structure samplelO. 
samplelO consists of all of the photometry for all of the targets over that area (as extracted from the internal 
SDSS operational database), all of the spectroscopic results (as output from specBS), and, most significantly, 
a description of the angular window function of the survey and the flux and surface brightness limits used for 
galaxies in each area of the sky. For most of the area, the same version of the image analysis software used 
to create the spectroscopic target list was used in this sample. However, for the area covered by the Early 
Data Release (EDR; Stoughton et al. 2002) we used the version of the analysis software used for that data 
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release, since it was substantially better than the early versions of the software used to target that area. For 
photo, the most important piece of analysis software run on the data, the versions used for the photometry 
range from v5_0 to v5_2. The region covered by this sample is similar to, but not exactly, the region which 
will be released in the SDSS Data Release 1 (DR1), scheduled for January 2003 (which will use photo v5_3, 
a newer version of the software which among other things improves the handling of large galaxies). Figure 
1 shows the distribution in right ascension and redshift of Main Sample galaxies with redshifts in samplelO 
within 6° of the Equator. Figure 2 shows the distribution of Main Sample galaxies with redshifts on the sky 
in Galactic coordinates. 

We measure galaxy magnitudes through a set of bandpasses which is obviously constant in the observer 
frame. This set of observer frame magnitudes corresponds to a different set of rest-frame magnitudes de- 
pending on the galaxy redshift. In order to compare galaxies observed at different redshifts, we convert all 
the magnitudes to a single fixed set of bandpasses. To perform this conversion, we use the method of Blanton 
et al. (2002a) (kcorrect vl_ll). These routines fit a linear combination of four spectral templates to each 
set of five magnitudes, assigning coefficients ao, a±, 02, and 04. The coefficient 00 to the first template is an 
estimate of the flux in the optical range (3500A < A < 7500A) in ergs s 1 cm 2 ; the fractional contribution 
of the other coefficients ai/ao, 02/00, and 03/00 characterize the spectral energy distribution of the galaxy. 
The most significant variation is along 03/ao. Taking the sum of the templates and projecting it onto filter 
responses, we can calculate the if -corrections from the observed bandpass to any rest-frame bandpass. In 
order to minimize the uncertainties in the if-correction, we choose a fixed set of bandpasses blucshifted by 
z = 0.1 in order that they cover the same region of the rest-frame spectrum as do the observed bandpasses 
for a galaxy at z = 0.1 (chosen because it is near the median redshift of the sample). The bottom panel 
Figure 3 shows this set of shifted bandpasses in the SDSS 01 g, °'V, °-H, 01 z)); the top panel shows the 

more commonly used unshifted restframe bandpasses, measurements of which are more poorly constrained 
by our set of data. 

For this paper, rather than using the full freedom of all four templates, we instead use if -corrections 
from a restricted set of models. We fix ai/o = and o 2 /o = (near the mean of the galaxy distribution 
in coefficient space). In addition, after fitting for 03/00 we restrict that ratio to be one of the twelve values. 
Figure 4 shows the if-corrections in each band as a function of redshift for each of these twelve values. The 
reason to restrict the ^-correction this way is to make the likelihood evaluation more efficient, as described 
in Blanton ct al. (2002c). We have tested the effect of this approximation on our fit for the luminosity 
function by taking a smaller set of values (six), finding no significant differences. Furthermore, if we consider 
the residuals of the restricted if-corrections in Figure 4 to the if-corrections found using all four templates, 
there is no trend with redshift and the standard deviations of the differences arc 0.03 mags at most (for the 
u band; the value is much less for the other bands). 

From this sample, we create one subsample for each band, u, <?, r, i, and z, of galaxies satisfying the 
apparent magnitude and redshift limits listed in Table 1. We chose the same apparent magnitude limits 
for the bands other than r as Blanton et al. (2001) chose; that paper chose the limits by imposing the 
requirement that < 2% of galaxies brighter than the flux limit in the given band are fainter than r* = 17.6. 
By defining a separate magnitude-limited sample in each band, we avoid biasing our results in one band 
due to the fact that the galaxies were selected in another band (for example, if we calculated the luminosity 
function in the u band from a sample limited in the r band, all of the faintest u band galaxies would tend 
to be red galaxies). For each band we include essentially all observed absolute magnitudes in our samples. 
However, we of course do not have good constraints at all absolute magnitudes; for this reason, we consider 
our model applicable only to a smaller range of evolution-corrected absolute magnitudes, which we define 
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from the tenth least luminous object to the most luminous object. These absolute magnitude limits are also 
listed in Table 1. 



3. Fitting the Luminosity Function 

In this section, we describe how to recover the number density of galaxies 3>(M, z) as a function of 
absolute magnitude M and redshift z, which we will use to calculate the luminosity density at z — 0.1. Each 
galaxy has a measured magnitude in each band m, an associated uncertainty Am, and a redshift z. The 
absolute magnitude Mo.i,, may be constructed from the apparent magnitude m r and redshift z as follows 

Mo.i r - 51og 10 h = m r - DM(z,ri ,fi A ,/i = 1) - Ka.i rr (z), (4) 

where as written DM(z,Qo,£Ia, h = 1) is the distance modulus as determined from the redshift assuming a 
particular cosmology (for example, using the formulae compiled by Hogg 1999) and h = 1. Ko.ir^z) is the 
if -correction from the r band of a galaxy at redshift z to the ° V band. 

In order to fit the distribution of the absolute magnitudes and redshifts of galaxies, we use a maximum 
likelihood method which allows for a generic luminosity function shape (it does not assume a Schechter 
function or any other simple form), allows for simple luminosity and number evolution, and accounts for the 
estimated uncertainties in the galaxy fluxes. The method is described in detail by Blanton et al. (2002c), and 
we outline it briefly here. It is akin to the stepwise maximum likelihood (SWML) method of Efstathiou et al. 
(1988). However, it includes evolution in our model for the luminosity function, it accounts for the effects of 
flux uncertainties, and (for computational convenience) it uses Gaussian basis functions rather than top hat 
basis functions. In addition, rather than maximizing the likelihood of absolute magnitude given redshift, it 
maximizes the joint likelihood of absolute magnitude and redshift. This choice makes our estimates more 
sensitive to large-scale structure in the sample, and more sensitive to evolution. 



Our model for the luminosity-redshift function is 

r 1 (M - M fc + (z - z )Q) 2 ] 



$(M, z) = nl0°- 4(z - zo)p V $ fc 1 exp 

k V 27rcr M 



a 2 



M 



where the Mk are fixed to be equally spaced in absolute magnitude and represent the centers of Gaussians 
of width om- &k are adjustable parameters signifying the amplitudes of the Gaussians. 

Following Lin et al. (1999), Q represents the evolution in luminosity, in units of magnitude per unit 
redshift; Q > indicates that galaxies were more luminous in the past. P quantifies the change in the number 
density with redshift; we choose this particular parametrization (again following Lin et al. 1999) such that P 
represents the contribution of number density evolution to the evolution in the luminosity density in units of 
magnitudes. P can be interpreted as either due to true evolution in the number density or due to very large 
scale structure. Given the size of our data set and its relative shallowness, we cannot distinguish between 
these possibilities; when necessary to, we will interpret P only as large-scale structure. In any case, our 
main interest in this paper is the luminosity density at z — 0.1, not its evolution, and the luminosity density 
is insensitive to reasonable values of P. zq is the zeropoint redshift, with respect to which we measure the 
evolution; for this sample, we choose zq = 0.1, the median redshift, since it is at that redshift that we can 
observe galaxies in the luminosity range around M*, which contribute the most to the luminosity density. 

In principle we can include large-scale structure in the radial direction, p(z), explicitly in the model, with 
the constraint that it have a reasonable power spectrum (since the power in the density field is constrained 
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mostly by modes which are not purely radial). However, we have decided not to do so here because it is not 
necessary for our goals. 

As described in Blanton et al. (2002c), we fit model parameters by maximizing the likelihood of the 
model parameters given the data: 



p(Q,P,ln<P k \M t ,z t ) = [[ pjM~7~) (6) 



We assume a uniform prior distribution of ln$fc, Q, and P (thus guaranteeing that is positive). Because 
p{Mi, Zi) obviously does not depend on the model parameters, the problem reduces to minimizing 

E= -2^1ogp(M J ,z J |Q,P,ln$ fe ). (7) 

We construct the likelihood p(M i; Zj \Q, P, In $&) of each galaxy i by convolving the luminosity function with 
a Gaussian of width Am (the estimated apparent magnitude uncertainty defined above) and constraining 
the galaxies to satisfy the flux limits of the survey: 

p(M z\Q P In $ fe ) = / $(M ' Z) ® G(Am) if mmin <M + DM(Z) + K(Z) ~ (Z ~ Zo)Q K mmax (8) 
' ' ' \ otherwise 

The number of parameters required for this fit (50-100) is small enough that standard function min- 
imizers can handle the task in a reasonable amount of time (one hour) on modern workstations (in our 
case, a 2 GHz Pentium IV machine), for a sample of ~ 10 5 objects. In the fit, we constrain the integral of 
<f>(M, z = 0.1) to be unity over our range of absolute magnitude (as listed in Table 1 for each band). 

The overall normalization n cannot be determined from this likelihood maximization procedure. We 
use the standard minimum variance estimator of Davis & Huchra (1982) to perform the normalization. 

n J dV^z)w(z) 1 [ > 

where the integral is over the volume covered by the survey between the minimum and maximum redshifts 
used for our estimate. The weight for each galaxy is 

W(z) = 7 -rS 7T, (10) 

V 7 l + nlO o - 4P (^ z «)J 3 0(z)' V ; 

and the selection function is 

<t>{z) = -, (11) 

J^dL^z) 

where ft is the galaxy sampling rate determined at each position of sky as the fraction of targets in each 
sector that were successfully assigned a classification. The integral of the correlation function is: 



J 3 = / drr 2 £(r) = 10000 Mpc 3 . (12) 

Clearly, because n appears in the weight w(z), it must be determined iteratively, which we do using the 
simple estimator n = (1/V) ^ ^/4>( z j) as an initial guess. 
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To determine the uncertainties in our fit, we use thirty jackknife resamplings of the data. In each 
sampling, we omit 1/30 of the effective area of sky (meaning, the area weighted by the sampling rate f t ). 
Each omitted area is a nearly contiguous set of sectors. This jackknife resampling procedure thus includes, 
to the extent possible, the uncertainties due to large-scale structure and calibration errors across the survey. 
Effectively, it includes the effects of errors which are correlated with angular position on the largest scales. 
Taking the results of all thirty fits to the data, we calculate the covariance between all of our measured 
parameters using the standard formula 

(AxiAxj) = — ^ J2( x * ~ Xi )( x i ~ x i) ( 13 ) 

The uncertainty correlation matrix is then defined in the standard way: = (AxiAxj) / ((Axf) (Ax: 2 )) 1 / 2 . 

We are not interested in the uncertainty correlations between the amplitudes of each Gaussian, because 
obviously neighboring Gaussians will be highly covariant. For this paper, we will not even be interested in 
the covariances between overall amplitude at different luminosities of the luminosity function calculated from 
the sum of the Gaussians. However, we will list in tables the covariances between the luminosity density, the 
evolution parameters, and overall measures of the shape, such as M* and a for the best fit Schechter function 
of each luminosity function. It is also important to track the covariances among the luminosity densities in 
all of the bands; because large-scale structure is an important source of uncertainty, the luminosity densities 
are highly covariant, and ignoring this covariance would lead to overconfidcnce in any fit to the stellar density 
in galaxies based on this data (and underconfidence in our knowledge of the relative luminosity density in 
different bands). 

Note that the true uncertainties in the luminosity density may be dominated by the uncertainties in 
the overall photometric calibration or by the fraction of flux contained within the Petrosian aperture for the 
galaxies that contribute to the luminosity density, while the uncertainties in the level of evolution recovered 
may be dominated by possible systematic errors in the if -corrections, as well as a systematic dependence of 
the fraction of light contained within a Petrosian magnitude as a function of redshift (due to the effects of 
seeing). 

4. Results 

We have applied the procedure described in the Section 3 to the data described in Section 2. Our results 
are summarized in this section. 



4.1. Luminosity Functions 

Figures 5 and 6 show the galaxy luminosity function in the 0A u, 01 g, °'V, °'H, and 01 z bands, assuming 
Oo = 0.3 and £1a = 0.7. The thick black line shows our best fit luminosity function. The thin black lines 
show the Gaussians which sum to form the full luminosity function. The grey region surrounding the thick 
black line indicates the la uncertainties in the luminosity function — of course, these uncertainties are all 
correlated with one another, and are closer to representing the uncertainties in the overall normalization of 
the function than the individual uncertainties at each magnitude. The best-fit Q and P evolution parameters 
are listed in the figure. 



-10- 



We have taken the thick black lines and their uncertainties and fit a Schechter function to each curve. 
The dotted lines in Figures 5 and 6 represent the best fit Schechter functions, which provide a reasonable 
fit to our non-parametric results. The luminosity density we list in the figure, expressed as the absolute 
magnitude from galaxies in an h~ 3 Mpc 3 on average, is the result of integrating this Schechter function fit 
over all luminosities. The values associated with the Schechter function are listed in Table 2. We list results 
for the (Clo = 1-0, &a = 0.0) and (fl = 0.3, Qa = 0.0) cosmologies as well. We have found that to an 
accuracy of about 3%, we can convert the results of one cosmology to those of another by scaling <p* by the 
inverse ratio of the comoving volumes at z = 0.1 between the two cosmologies, and by scaling M* by the 
difference of the distance moduli at z = 0.1 for the two cosmologies. We therefore recommend this procedure 
for readers interested in comparing our results to those in some other cosmological model. 

Table 3 lists some salient quantitative measurements of the luminosity function in each band, including 
the evolution parameters and the luminosity density (expressed in magnitudes, solar luminosities, and flux at 
the effective filter wavelength) for a Mpc 3 . To obtain the physical expressions of the luminosity density, we 
used measurements of the SDSS camera response performed by Mamoru Doi, which James Gunn combined 
with estimates of the atmospheric extinction as a function of wavelength at 1.3 airmasses (to which all SDSS 
observations are calibrated) and the primary and secondary mirror reflectivities. Projecting the solar model 
of Kurucz (1991) onto these bandpasses (shifted to z = 0.1) yields the absolute solar AB magnitudes: 

M Q o.i, a = 6.80; M 0i o.i g = 5.45; M 0> o. v = 4.76: M e o. H = 4.58; M 0; o.i 2 = 4.51 (14) 

The luminosity densities expressed in ergs s _1 A -1 are calculated from the AB magnitudes as follows. First, 
we use the equation which relates an AB magnitude to the effective flux density at the effective wavelength, 

fx = (3.631 x 10 - 20 ergs cm" 2 s" 1 Hz" 1 )-^^ - 4 ™ (15) 

Ks 

to convert the absolute magnitudes in an h~ s Mpc 3 to the flux density which would be observed if an average 
h~ 3 Mpc 3 of the universe were compressed to a point source and placed 10 pc distant from the observer. 
Second, we multiply this value by 47r(10 pc) 2 to obtain the average luminosity per unit wavelength at the 
effective wavelength in an h~ 3 Mpc 3 . The "effective wavelength" of a passband with a quantum efficiency 
R(X) is defined: 

"/ d{\nX)R{X) In A" 



A e ff = exp 



(16) 



/ d(ln\)R(X) 

following Fukugita et al. (1996) and Schneider et al. (1983). The effective flux density defined above is that 
which an AB standard source (fx (A) cx A~ 2 ) of magnitude m in passband R would have at the effective 
wavelength. Both of these quantities are obviously only nominal since, in any case, the average spectrum of 
galaxies is nothing like an AB standard source, but it does give a sense of the physical flux associated with 
a magnitude. We also list / np , the fraction of the integrated luminosity density of the Schechter luminosity 
function included in the non-parametric estimate of the luminosity density. Results for the cosmologies 
(Cl m = 0.3, Oa = 0.0) and (Q TO = 1.0, Oa = 0.0) are also listed. 

It is worth asking how well this model reproduces the number counts of galaxies as a function of rcdshift 
and absolute magnitude. Figure 7 shows the redshift distribution of galaxies in our sample for quartiles 
in absolute magnitude in the ° V band as a thick histogram. The expectation from our model fit (based 
on Monte Carlo realizations of the sample) is shown as the thin histogram. Figures 8-12 show the counts 
in bins of absolute magnitude for several slices in redshift. Again, the thick histogram represents the data 
and the thin histogram represents the model. The model appears to the eye to be a reasonable fit to the 
data. However, it is clear that there are statistically significant discrepancies in these figures (note the large 
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number of objects in each bin). Some of these discrepancies occur because there is large-scale structure 
in the sample; however, it is possible that we will need to introduce a more sophisticated model for the 
evolution of the galaxies or for the relationship between density and luminosity in order to explain all of 
these discrepancies. As discussed in the conclusions, we postpone the investigation of these issues to a later 
paper. 

As described above, we have quantified our uncertainties by taking 30 jacknife resamplings of the data. 
In Tables 4-8 we display the resulting correlation matrices between various properties of our fit for each 
band. Note that many of the parameters are highly correlated. In particular we note that M* and a, which 
characterize the shape of the Schechter luminosity function, are highly correlated. One should be cautious 
when quoting M* values as representative of the "typical luminosity" of a luminous galaxy. This statement 
is true to an extent, but not in detail. 

Table 9 displays the correlation matrix between the luminosity densities and evolution parameters Q in 
all of the bands. Note that this matrix is entirely positive. In any particular resampling, if the luminosity 
density is slightly high in one band relative to the mean, all the bands are slightly high; when one band 
evolves a bit more strongly, so do they all. This correlation occurs because our errors are dominated by large- 
scale structure, which affects all the bands simultaneously. These correlations are strong, so it is important 
to account for this correlation matrix of uncertainties when using these results. 

The top panel of Figure 14 shows the luminosity density as a function of wavelength at z = 0.1 using 
the results of all five bands. 



4.2. Galaxy Luminosity Density Evolution 

One product of our fit is Q, the evolution of the luminosity density. However, we caution that we have 
restricted our model such that galaxies of all luminosities evolve identically. This assumption is probably 
incorrect, because different galaxy types have different luminosity function and are expected to evolve in 
different ways (since their differing colors obviously imply different star- formation histories). Furthermore, 
our understanding of our photometric error model, particularly in u and z, is currently rather primitive. 
While we believe that the flaws in our model for the evolution do not greatly bias the main result of this 
paper, the luminosity function and luminosity density at z = 0.1, we do warn the reader that because of the 
deficiencies of our model the measured Q values might be biased. 

We show our results for the luminosity evolution, Q, and its uncertainties, in the bottom panel of Figure 
14, for all the bands. The band has very strong evolution (Q ~ 4), although with large uncertainties. 
The other bands all have Q ~ 1-2. The evolution in these bands is generally consistent with the evolution 
of a relatively old stellar population. However, we will consider fits to theoretical models more carefully in 
a separate paper, and draw no particular conclusions here. 

It is likely that these evolution results reflect the evolution of the more luminous galaxies in the sample, 
for two reasons. First, more luminous galaxies are observable over a larger redshift range. Second, the low 
luminosity end of the luminosity function is nearly a power-law, which makes it more difficult to detect 
evolution. If Q traces the evolution of primarily the luminous galaxies, it might explain why the evolution in 
the band is unusually large; the most luminous objects in the ° \t band are the blue, exponential profile 
objects, which we expect to evolve more rapidly than red, concentrated galaxies. 

For the flat, A-dominated cosmology, o l u and 01 z bands have P^O, though at low significance in the 
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case of 01 u. We do not believe that P reflects true number density evolution, which we regard as a priori 
unlikely, especially since the other bands have P consistent with zero. Instead P probably is compensating 
for inaccuracies in some other aspect of our model. In the case of the o l u band it may compensate for the 
fact that our model does not include large-scale structure. In the case of the 01 z band, our error model 
could be insufficient, or the if-corrections might be incorrect (though they are small in any case). Thus, 
one might claim that the luminosity density evolution we measure should be Q + P; however, we think it 
better to interpret a non-zero P to reflect large-scale structure in the sample rather than true evolution. A 
more sophisticated approach would be to treat the problem in more detail by fitting for the radial large-scale 
structure simultaneously with the luminosity function. 

We note in passing that the differences in P and Q between the flat, A-dominated cosmology and the 
other cosmologies are consistent with the differences in the redshift dependence of the luminosity distance 
and the differential comoving volume among these cosmologies. We do not regard the low value of P for g, 
r, and i in the flat, A-dominated cosmology necessarily as evidence that that cosmological model is correct. 

5. Systematics Tests of Petrosian Magnitudes 

One of the challenges in studying the luminosity function and its redshift dependence is that one must 
verify that one's measurements of galaxy luminosity are consistent as a function of redshift. That is, one 
must test whether observing an identical galaxy at different rcdshifts will yield the same fraction of the 
total independent of redshift. Otherwise, systematic errors in the measurement of galaxy luminosities as a 
function of redshift could masquerade as evolution. 

Several :w effects can yield artificial redshift dependence in one's determination of galaxy flux (which is 
to say, redshift dependence which cannot be accounted for using the distance modulus and the if -correction) . 
First, the increased physical size corresponding to the angular PSF width for galaxies at high redshift can 
change both the metric size (in physical units at the galaxy) of the aperture used to calculate the galaxy flux 
and the amount of light scattered outside of any particular aperture. Second, color gradients in galaxies mean 
that the different radii in the galaxy profile if-correct differently. Since in the SDSS the size of the aperture 
is based on the shape of the profile, this effect can change the metric size of the aperture used to calculate the 
galaxy flux. Third, cosmological surface brightness dimming can make the outermost measurable isophote 
have a smaller metric size. 

Because the SDSS has high signal-to-noise ratio imaging and because it uses Petrosian rather than 
isophotal magnitudes, the third effect is negligible. However, Petrosian magnitudes are sensitive to some 
degree to seeing; as an object becomes small with respect to the seeing width, the fraction of the total flux 
measured by the Petrosian magnitudes tends to that fraction for a pure PSF (e.g. Blanton et al. 2001, 
Strauss et al. 2002). For the SDSS, the fraction of the total light in a pure PSF measured by the Petrosian 
magnitude is about 95%. Thus, galaxies with exponential profiles will tend to have their measured fluxes 
reduced in the presence of bad seeing from the fraction of nearly 100% for the no-seeing case. On the 
other hand, galaxies with de Vaucouleurs profiles will tend to have their measured fluxes increased from the 
fraction of only 82% for the no seeing case. Furthermore, the apertures of the Petrosian magnitudes are 
determined by the r-band image of the galaxy, which, while it is less sensitive to differential if-corrections 
than bluer bands, is still somewhat dependent on these effects. 

Several methods exist for estimating what fraction of the flux is outside the SDSS Petrosian apertures. 
First, one could compare SDSS imaging to deeper imaging. Second, one could stack many SDSS images 
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of different galaxies but of similar types, and use the composite profile to characterize the fraction of light 
missing outside of a Petrosian aperture. Third, one can use a reasonable model for galaxy light profiles to 
extrapolate to a "total" flux for each object. Here, we take the third approach, using the Sersic magnitudes 
described in Section 2 as a seeing-corrected, "total" magnitude. In addition, because the Sersic fit is based 
on the i-band image, the issue of the Zf-correction of the galaxy profile is negligible (as the if-corrections in 
the z-band are small) for the redshifts considered in this paper. 

Figure 13 shows the differences between Sersic and Petrosian magnitudes in each band, as a function of 
Sersic index n and redshift z. Clearly, the difference is sensitive to the Sersic index and is larger for galaxies 
which are similar to de Vaucouleurs galaxies, in accordance with the estimates in Blanton ct al. (2001) of 
the fraction of flux included in the SDSS Petrosian magnitudes for different profile types. For u and g the 
Sersic magnitude is significantly fainter than the Petrosian magnitude at low Sersic index. This fact might 
reflect color gradients within exponential galaxies. The shape of the Sersic model we use is determined in the 
i-band; if the effective size of the galaxy is larger in bluer bands, known to be the case for spiral galaxies, the 
flux determined by fitting the Sersic model amplitude will be smaller than the flux determined by counting 
the flux within an unweighted aperture. The right panel shows a linear regression versus redshift. The mean 
offset at z = 0.1 ranges from 0.14 for the ° hi band to -0.14 for the 0A z band. Again, this result suggests that 
the Sersic magnitude traces a redder population, perhaps as a consequence of color gradients in galaxies. 

For our purposes here, we are interested primarily in the differences in the resulting luminosity density 
from Sersic magnitudes relative to Petrosian magnitudes. When we fit luminosity functions to all five bands 
using the Sersic magnitudes, we obtain luminosity densities offset from the Petrosian values as follows: 
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The values of Q for Sersic magnitudes are consistent with the results for Petrosian magnitudes. We caution 
that the estimates of the luminosity density using Sersic magnitudes are not necessarily more correct than 
the estimates from Petrosian magnitudes; because the Sersic profile is not a perfect fit to the observed radial 
profiles of galaxies, there is likely to be a bias in the resulting estimates of luminosity density associated 
with Sersic magnitudes. However, the small differences in Equation (17) suggest that the differences are 
also small between either of these estimates and the luminosity density we would derive using true "total" 
magnitudes. 

6. Comparison with Other Results 

6.1. Galaxy Luminosity Density in Other Bands 

In order to compare our results to those of other investigators, we have used the routines in kcorrect 
vl_14 (a slight update on the version kcorrect vl_ll used elsewhere in this paper and in Blanton et al. 
2002a) to fit the SED of the universe and derive the luminosity densities in other bands. In particular, we 
list in Table 10 results for the unshifted SDSS bandpasses for comparison to Blanton ct al. (2001), for the 
unshiftcd B, V, R and / Johnson bandpasses (using the responses listed by Bessell 1990), and for bj for 
comparison to 2dFGRS. To create the bj result we applied the equation bj—B = —0.28(5 — V) to the results 
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in the Johnson bandpasses, originally from Blair & Gilmore (1982). We derive these luminosity densities 
both at z = 0.1 and (by using Q to evolve the SDSS results) z = 0. For comparison, we show in the same 
table the results of Blanton et al. (2001) and Norberg et al. (2002). For the Norberg et al. (2002) we use 
the mean evolutionary correction from their Figure 8 to evaluate a luminosity density at z = 0.1. There is 
only 0.09 magnitude difference between our result in bj at z = 0.1 and theirs; this difference is less than la 
taking into account the 2dFGRS statistical uncertainties. 

The results of Table 10 reflect a general agreement between these different determinations of the lumi- 
nosity function. In contrast to these results, Blanton et al. (2001) reported a higher luminosity density than 
found in other surveys, such as the 2dFGRS results of Folkes et al. (1999) and the Las Campanas Redshift 
Survey (LCRS; Shectman et al. 1996) results of Lin et al. (1996). Most recently, Liske et al. (2002) found 
that the 00 <7-band luminosity function of Blanton et al. (2001) overpredicted the number counts of galaxies 
found in the Millenium Galaxy Catalog. Furthermore, as pointed out by Wright (2001), the luminosity 
density in the z-band reported by Blanton et al. (2001) was unreasonably high relative to the luminosity 
density found by Cole et al. (2000) in J and 00 K. While our luminosity density is still considerably higher 
than that found by the LCRS, we no longer disagree significantly with 2dFGRS. Naturally, the question 
arises as to why the discrepancy existed in the first place. 

Table 10 makes clear that there are large differences between our present results and those of Blanton 
et al. (2001). We should emphasize here that if we apply the same methods used in that paper to our current 
sample, the results are within la of the results in that paper. This change in the luminosity density estimates 
are not due to statistical fluctuations, nor to a change in the photometric calibration, but are instead due to 
improvements in our model for the luminosity function and the observations. In particular, we are using a 
better estimate of the -corrections and we are including evolution in our model for the luminosity function. 
The largest changes are the luminosity densities in the u and g bands (which are 1.2 and 0.75 magnitudes 
less luminous respectively), where the changes in the X-corrections were largest and where evolution is most 
important. Ignoring evolution in the luminosity function model causes a large bias in the estimate of the 
luminosity density, because it causes the expected number of objects at high redshift to be inaccurate. If 
galaxies in fact are brighter in the past, a non-evolving model tends to yield lower number counts at high 
redshift, at a given normalization. Since the normalization procedure of Davis & Huchra (1982) weights 
according to volume, and thus accords higher weight at higher redshift, in this case a non-evolving model 
would result in an overestimate of galaxies at low redshift. Due to bad luck, the systematics comparison 
of Figure 8 in that paper, which compared the normalizations of the luminosity function at high and low 
redshift, happened not to reveal this effect, presumably due to the large supercluster at z ~ 0.08 in those 
data (and still distinctly visible in Figure 7 in this, much larger, dataset!). Figures 7-12 in the present paper 
show decisively that our current model explains the redshift counts very well. 

So how does this affect our comparisons to other surveys? For the LCRS, whose method of fitting 
the luminosity function and its normalization was identical to that of Blanton et al. (2001), the original 
comparison remains the fair one. That is, even though our estimate of the luminosity density is now only 0.2 
magnitudes brighter than that of Lin et al. (1996), this is only an accident, resulting from a combination of 
two effects in the LCRS: using bright isophotal magnitudes, which lowers the luminosity density estimate, 
and ignoring (as Blanton et al. 2001 also did) evolution, which raises the luminosity density estimate. For 
2MASS, the change in our result makes the SDSS more compatible with the results of Cole et al. (2001). 
However, it is more difficult to directly compare these surveys, since the SDSS bands and 2MASS bands do 
not overlap. 

For the 2dFGRS, Norberg et al. (2002) report a luminosity density of jb } = —15.35 absolute magnitudes 
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at z = (integrating the Schechter function for the fio = 0.3, fl\ = 0.7 cosmology in the first line of their 
Table 2 over all luminosities). This result is based on extrapolating to z — the luminosities of galaxies 
whose typical redshifts are z = 0.05-0.2, using assumptions about the luminosity evolution. Figure 8 of 
Norberg et al. (2002), which shows the mean if-correction and evolution correction used in their analysis, 
shows that their evolution correction corresponds closely to Q = 1. Since we find a somewhat different 
value of Q ~ 2 at these wavelengths, and both surveys have similar median redshifts, the fair comparison 
of the luminosity densities involves evaluating the luminosity density at around z = 0.1. For this reason, 
we evolution-correct their results back to z = 0.1 by applying Am = — 0.1Q = —0.1. Thus, for 2dFGRS 
jbj (z = 0.1) = —15.45 ± 0.1, within la of our result in Table 10. Note that if we instead compare our z = 
value of the bj luminosity density to theirs, the discrepancy is about 0.2 magnitudes. However, in cither 
comparison the differences between the SDSS and 2dFGRS luminosity densities are rather small. 

We find similarly small differences between our results and the galaxy counts data from the Millenium 
Galaxy Catalog of Liske et al. (2002). Our Schechter function model for the al g-band luminosity function, 
crudely transformed into -Bmgc bandpass by applying —0.09 to M», and using Q = 2 as for 01 g, nicely 
predicts their galaxy counts to -Bmgc — 20 to within about 5%. 

The original Folkes et al. (1999) 2dFGRS result that Blanton et al. (2001) compared to did not evolution- 
correct their magnitudes either. So why was there a discrepancy between those two results? First of all, as 
Norberg et al. (2002) have pointed out, the linear transformation used by Blanton ct al. (2001) to convert 
SDSS magnitudes to 2dFGRS magnitudes (which was based on the linear transformation between bj and B 
and V published by Metcalfe et al. 1995) was inappropriate. Using the Blair & Gilmore (1982) transformation 
from B and V to bj (bj —B = —0.28(5 — V)) instead results in a 10% reduction in the luminosity density, not 
enough to explain the original discrepancy. However, one significant difference between Folkes et al. (1999) 
and the mock bj luminosity function in Blanton et al. (2001) was that Folkes et al. (1999) normalized to 
the number counts instead of using a volume-weighted method. Due to this normalization, the Folkes ct al. 
(1999) result was less affected by luminosity evolution, because their normalization was set by galaxies close 
to the median redshift of the survey, rather than those at the distant edge of the survey. This difference 
seems to account for the bulk of the discrepancy between Blanton et al. (2001) and Folkes et al. (1999). 
Norberg et al. (2002) reach the same conclusion. As stated above, comparison of our results to those of 
Norberg et al. (2002), who perform an evolution correction to their magnitudes, now results in a consistent 
measurement of the luminosity density. 

Nevertheless, Blanton ct al. (2001) did show that the isophotal magnitudes meaured by the APM 
survey did not include a significant fraction of the flux for a typical galaxy, and further that the "Gaussian- 
correction" to total magnitudes (described in detail by Maddox et al. 1990) performed on the galaxy fluxes 
were insufficient to replace this difference. If this were the whole story, there should remain a difference 
of roughly 30% between the Norberg et al. (2002) results based on 2dFGRS magnitudes and the results 
presented here, based on SDSS Petrosian magnitudes. However, extensive further calibration of the plate 
magnitudes has been performed using deep CCD imaging, as described by Norberg et al. (2002). In effect, 
the overall calibration of the galaxy survey is set by a comparison between observations of a set of galaxies 
using the APM and using deep CCD imaging. Since the magnitudes in the deep CCD imaging count flux 
outside the limiting isophotes of the APM magnitudes, this overall calibration translates the isophotal to 
total magnitudes (on average, that is — naturally, for any particular galaxy the correction depends on the 
profile shape and size of that galaxy). For this reason, the comparisons in Norberg et al. (2002) of 2dFGRS 
magnitudes and SDSS magnitudes agree on average, for galaxies near the median surface brightness (although 
there is a strong correlation of the residuals with surface brightness). Blanton et al. (2001) did not account 
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for these extra step of calibration and isophotal-to-total correction, and so overestimated the effect of the 
different magnitude definitions in the final luminosity density. 

Both Liske et al. (2002) and Cole et al. (2001) have suggested that the results of Blanton et al. (2001) 
could have been biased high due to large-scale structure. This possibility can be ruled out by simply using 
the same methods used by that paper on the full sample presented here. The luminosity density so calculated 
is within lcr of the original result, suggesting that the 10% uncertainties calculated for that sample were 
realistic. This agreement occurs despite the fact that the galaxy counts for r < 17.7 in the region used 
by Blanton et al. (2001) are 10-15% higher than in the survey on average. As correctly pointed out by 
Norberg et al. (2002), the statistical uncertainties in the normalization arc much smaller when using the 
volume-weighting method of Davis & Huchra (1982), rather than the normalizing to galaxy counts, though 
as those authors state one is more susceptible to systematic errors in the model for evolution. Using the 
luminosity function model in this paper, normalizing either way yields nearly identical results; for example, 
consider the counts predicted by the models listed in Figures 8-12. The possible exception is the w-band, for 
which there is a 10% difference (still within the la statistical uncertainties). Again, we emphasize that the 
differences between the luminosity function derived here and in Blanton et al. (2001) are a product primarily 
of the model of the luminosity function we use. 

In short, the discrepancies between SDSS and 2dFGRS appear to be resolved, for the most part due to 
a proper treatment of evolution in the SDSS luminosity function. It remains to be determined whether the 
results of SDSS and 2MASS are discrepant, but in any case they are ameliorated by the reduction of our 
estimate of the luminosity density. 

6.2. Luminosity Density as a Function of Wavelength 

Figure 14 plots the luminosity density versus wavelength using Equation (15) for the SDSS (diamonds). 
Using the above results, we add the 2dFGRS point from Norberg et al. (2002) (triangle), and the 2MASS 
points from Cole et al. (2001) (crosses), as described in this section. 

For the 2dFGRS, translating the luminosity density expressed in magnitudes into physical units is a 
nontrivial endeavour. The zcropoint of the 2dFGRS bj system is set by the equation: 

bj — B = — 0.28(S — V). (18) 

Although this slope was determined by Blair & Gilmore (1982) based only on a small set of stars, it can 
be checked a posteriori against CCD observations which are appropriately converted from their natural 
system into the Johnson B and V system. According to Norberg et al. (2002), such observations confirm 
the slope —0.28 to a surprising degree of accuracy given the estimated errors of Blair & Gilmore (1982). 
However, the typical galaxy used to perform this check is not anywhere near B — V = 0, and because we 
expect the relationship between any set of bandpasses to be nonlinear (with or without photographic plate 
nonlinearities), it is not guaranteed that bj, as zeropointed using Equation (18, is a Vega-relative magnitude. 
We perform a test of Equation (18) using an estimate of the spectral response of the plates, from Hcwett & 
Warren (private communication), who have determined the response for one plate by dispersing a spectrum 
(with well- measured f\ (A) from spectrophotometry) through 2mm of GG395 filter and onto Kodak IllaJ 
emulsion, the definition of the bj bandpass. Using the nominal B and V responses derived by Bessell (1990) 
and typical galaxy spectra at z — 0, we find that the relationship between bj and B around a typical galaxy 
color of B — V <~ 0.8 is close to bj—B « 0.04 — 0.28(i3 — V) . The relationship is indeed nonlinear and appears 
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to cross bj — B = at B — V = 0. Again, if this difference existed in the 2dFGRS data, it would never be 
noticed because the zeropoint is explicitly set using CCD imaging of galaxies assuming Equation (18) holds 
at typical galaxy colors. On the other hand, the curve given to us by Hewett & Warren is based only on 
measurements of a single plate, whose response may not be typical of the plates in the APM. Given this 
uncertainty, we simply assume Hewett & Warren are correct, and derive an AB magnitude from 2dFGRS bj 
magnitudes: 

bj (AB) = bj (2dFGRS) + 0.04 - 0.07, (19) 

which we claim is uncertain at the 0.05 mag level. Then we can use Equation (15) to calculate the luminosity 
density in physical units. 18 

In the near infrared, the Two-Micron All Sky Survey (2MASS) results of Cole et al. (2000) for the 
°-°J and K S bands are shown. We use the published responses of J and K 8 from the 2MASS website 19 
and the Kurucz (1991) theoretical Vega spectrum (normalized at 5000 A to match the Hayes (1985) spec- 
trophotometry of Vega) to calculate zeropoint shifts from Vega to AB magnitudes of (0.91, 1.85) for J and 
K s , respectively. The J band AB zeropoint corresponds closely with that listed in the 2MASS Explanatory 
Supplement Section IV. 5. a. However, the K s band AB zeropoint listed there is about 1.77; since the spec- 
trophotometry or model used to calculate this zeropoint is not specified on that site, we cannot determine 
the source of this discrepancy. However, the difference is within the la statistical uncertainties of Cole 
ct al. (2001). To evolution-correct the 2MASS luminosity densities to z = 0.1, we assume Q = 1, which is 
consistent with most stellar population synthesis models. We display this in Figure 14 as crosses. 

We should note here that although we have done our best here to place all the observations on the same 
physical footing, absolute calibration of this sort is uncertain. No spectrophotometry we know of has been 
performed to verify the models of Kurucz (1991) in the infrared, which we rely on to put the 2MASS results 
in physical units in Figure 14. There are uncertainties of at least 5% percent in the spectrophotometric 
calibration of BD+17 4708, the primary standard used to calibrate the SDSS. In any of the observations 
there is uncertainty and quite possibly variability in the bandpass responses. So this plot (and any plot like 
it in the literature) should not be trusted to better than 5% (which can be enough to substantially change 
one's theoretical interpretation of the observations in terms of a star- formation history). 

We leave until a later paper or to other authors a serious attempt to reconcile these various observations 
given a star-formation history and stellar population synthesis models. 



7. Conclusions and Future Work 

We have presented the galaxy luminosity density across the optical range, and have shown that it is 
consistent with other determinations. Previous differences found between the SDSS and other surveys are 
primarily attributable to the effect of luminosity evolution on our results. We have accounted for luminosity 
evolution in our current fits, and have given our results, which appear loosely consistent with the predictions 
of passive evolution in most bands. Our results are similar to those found by Bernardi et al. (2002) for their 
sample of early-type galaxies. 

We note that our results are in general agreement with the shape of the cosmic spectrum determined 



18 This procedure yields a conversion between magnitudes and physical units substantially different (by about 50%) than what 
one would infer about the conversion from Table 1 of Folkes et al. (1999). 

19 http: //www. ipac. caltech. edu/2mass/releases/second/doc/sec3_lbl .tbll2.html 
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from the 2dFGRS and SDSS galaxy spectra by Baldry et al. (2002) and Glazebrook et al., in preparation. 
In the case of SDSS spectra, Glazebrook et al. use the same data as used here (Petrosian magnitudes in 
samplelO) to apply an overall normalization to their spectrum, but use the 3" diameter fiber spectra (which 
typically contain about 20% of the Petrosian galaxy flux and cover regions of the galaxy about 0.1 mag 
redder in g — r than that covered by the Petrosian aperture) to find the luminosity density as a function of 
wavelength. The power of their approach is that the emission and absorption lines contain more detailed 
information on the star-formation history than do broad-band colors; the question remains whether this 
star- format ion history is representative of the global star-formation history. 

We will continue to improve our estimates of evolution. First, we will quantify our uncertainties in 
flux measurements more robustly and understand how our estimates of magnitude are affected by seeing in 
more detail. Doing so will yield more reliable estimates of luminosity evolution. Second, different galaxy 
types are expected to evolve differently; redder galaxies are generally presumed to be older and thus more 
slowly evolving. The different evolution of different galaxy types can constrain theories of the formation of 
galaxies, and thus it is of interest to study this differential evolution. Furthermore, it is well known that 
different galaxy types have different luminosity functions (e.g., Blanton et al. 2001); thus, it is probable that 
the shape of the galaxy luminosity function evolves, not merely its overall luminosity scale. 

One caveat to our results is the existence of surface-brightness selection effects. Blanton et al. (2001) 
demonstrated that if the SDSS can truly detect most of the objects which exist down to fj, r ^ ~ 24, not 
much luminosity density (~ 5%) can exist at lower surface brightnesses, unless the dependence of luminosity 
on surface brightness flattens dramatically below these limits or there is a sharp upturn in the luminosity 
function at low luminosities. Cross et al. (2001) demonstrated similar results for 2dFGRS. This result is not 
necessarily inconsistent with the results of O'Neil & Bothun (2000), who found that the number density of 
galaxies does not decline at low surface brightnesses; however, if the relationship between luminosity and 
surface brightness measured in the SDSS exists in their galaxy sample, the contribution to the luminosity 
density of galaxies should decline considerably at low surface brightness. Unfortunately, that paper and its 
predecessors do not evaluate the luminosity density contributed by different ranges of surface brightness in 
their sample. 

The possibility remains that the SDSS is very incomplete at fx r ^o ~ 24. Such incompleteness would 
probably not be due to low signal-to-noise (see the order-of-magnitude calculation in Blanton et al. 2001); it 
would more likely be due to systematic difficulties in subtracting the sky background. We find no evidence for 
changes in the measured extinction-corrected surface-brightness distribution at Galactic extinctions varying 
from to nearly 2 in the r-band (as determined by Schlcgcl et al. 1998), which indicates that the surface 
brightnesses of most galaxies in the survey are not close to the surface brightness limit of the survey. In fact, 
there are a number of "objects" detected at low surface brightness /i r ,5o > 24, nearly all of which turn out 
to be scattered light or other image defects; thus, such low surface brightness features are readily detectable. 
We can test these effects more thoroughly by searching for simulated galaxies inserted into actual data, and 
such tests are currently ongoing. 

These results (and the accompanying covariance matrices) can be used to constrain aspects of the star- 
formation history of the universe, and in particular the overall stellar density. In addition, they can be used 
to develop a selection function for flux-limited galaxy redshift surveys selected from the SDSS. They also 
provide the state-of-the-art luminosity densities with which to calculate TO using measured mass-to-light 
ratios. Finally, they can be compared to high redshift estimates of luminosity density in order to evaluate 
the evolution of galaxies to high redshift. 
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Table 1. Selection limits 



Band 


Flux Limit 


Rcdshift Limits 


Absolute Magnitude Limits 


Number of Galaxies 


°-hi 


m„ 


< 18.36 


0.02 < z < 0.14 


-21.93 < Mo.i^ 


- 5 log 10 h < 


-15.54 


22020 


0.1 g 


m g 


< 17.69 


0.02 < z < 0.17 


-23.38 < Ma.i g 


- 5 login h < 


-16.10 


53999 




m r 


< 17.79 


0.02 < z < 0.22 


-24.26 < Mo. v 


- 5 log 10 h < 


-16.11 


147986 


o. H 


m l 


< 16.91 


0.02 < z < 0.22 


-23.84 < Mo. H 


- 5 log 10 h < 


-17.07 


88239 


o.i z 


m z 


< 16.50 


0.02 < z < 0.22 


-24.08 < Mo.i z 


- 5 login h < 


-17.34 


73463 



Note. — Arbitrarily apparently bright objects were included as long as they passed the galaxy target 
selection limits of Strauss ct al. (2002). Absolute magnitude limits are those used for the (r2 = 0.3, 
£1a = 0.7) case. 



Table 2. Schechter function fits 



n 




Band 


0* (10- 2 /i 3 Mpc" 3 ) 


M* - 5 log 10 h 


a 


0.3 


0.7 




3.05 ±0.33 


-17.93 ±0.03 


-0.92 ±0.07 






0.1 g 


2.18 ±0.08 


-19.39 ±0.02 


-0.89 ±0.03 








1.49 ±0.04 


-20.44 ±0.01 


-1.05 ±0.01 








1.47 ±0.04 


-20.82 ±0.02 


-1.00 ±0.02 






o.i z 


1.35 ±0.04 


-21.18 ±0.02 


-1.08 ±0.02 


0.3 


0.0 




3.26 ±0.40 


-17.89 ±0.04 


-0.94 ±0.09 






0.1 g 


2.42 ±0.10 


-19.34 ±0.02 


-0.92 ±0.04 








1.69 ±0.06 


-20.37 ±0.02 


-1.03 ±0.03 








1.62 ±0.06 


-20.78 ±0.03 


-1.02 ±0.04 






o.i z 


1.47 ±0.05 


-21.12 ±0.02 


-1.07±0.03 


1.0 


0.0 




3.65 ±0.40 


-17.83 ±0.04 


-0.90 ±0.06 






0.1 g 


2.62 ±0.10 


-19.30 ±0.02 


-0.91 ±0.03 








1.83 ±0.06 


-20.33 ±0.03 


-1.04 ±0.03 








1.73 ±0.06 


-20.74 ±0.02 


-1.03 ±0.03 






o.i z 


1.57 ±0.05 


-21.11 ±0.02 


-1.10±0.02 



Note. — The uncertainties are correlated; see Tables 4-8 for the corre- 
lation matrix of the uncertainties. 



Table 3. Luminosity density and evolution parameters 







Band 


Aoff 
(A) 


j + 2.51og 10 /i 
(mags in Mpc 3 ) 


3 

(h 10 8 L Q Mpc" 3 ) 


3 

(h 10 37 ergs cm^Mpc" 3 ) 


Q 


P 




/np 


a 9 


U. ( 


O.L, 


ozlb 


— 14. 1U ± U.lo 


o on _i_ n oo 
z.zy ± U.oz 


0.4» ± U. (0 


A OO _L A QQ 

4.ZZ ± U.oo 


o on _1_ O 


O 1 

,ol 


A OA 

u.yu 






0.1„ 

9 


A O A A 

4z40 


1 C 1 o I n no 

— 15.18 ± 0.03 


1.78 ± 0.05 


o.54 ± O.zo 


o a a _i_ a c 1 
z.U4 ± 0.51 


a oo _i_ 1 
0.3z ± 1 


TA 

.70 


A AT 
0.97 








5595 


i c aa _l a no 
— 15.90 ± 0.03 


1 o a _i_ a r\ /( 
1.84 ± 0.04 


9.57 ± 0.22 


i £?o _i_ a on 
1.62 ± 0.30 


A 1 o _L a 

0.18 ± 


.57 


1 A A 
1.00 








6792 


-16.24 ±0.03 


2.12 ±0.05 


8.82 ±0.21 


1.61 ±0.43 


0.58 ±1 


.06 


0.99 






o.i z 


8111 


-16.56 ±0.02 


2.69 ±0.05 


8.33 ±0.15 


0.76 ±0.29 


2.28 ±0. 


.79 


1.01 


0.3 


0.0 




3216 


-14.14 ±0.19 


2.39 ±0.42 


5.70 ± 1.00 


3.67 ±0.89 


4.02 ±3 


.18 


0.90 






0.1 g 


4240 


-15.26 ±0.07 


1.92 ±0.12 


9.17 ±0.55 


1.22 ±0.73 


1.82 ±2 


.10 


0.97 






o.if. 


5595 


-15.95 ±0.03 


1.93 ±0.06 


10.03 ±0.30 


1.11 ±0.48 


0.99 ±0 


.96 


1.00 






°-H 


6792 


-16.31 ±0.04 


2.27 ±0.09 


9.47 ±0.36 


0.94 ± 0.46 


1.71 ± 1 


.14 


0.98 






o.i z 


8111 


-16.59 ±0.04 


2.75 ±0.09 


8.52 ±0.29 


0.48 ± 0.48 


2.54 ± 1 


.23 


1.01 


1.0 


0.0 




3216 


-14.18 ±0.14 


2.48 ±0.33 


5.94 ±0.78 


3.33 ±0.66 


4.89 ±2. 


.67 


0.90 






0.l g 


4240 


-15.29 ±0.05 


1.99 ±0.09 


9.50 ±0.42 


0.95 ±0.49 


2.52 ±1. 


.61 


0.97 








5595 


-16.01 ±0.04 


2.03 ±0.07 


10.54 ±0.34 


0.63 ±0.55 


2.16 ±1. 


.07 


1.00 






°-H 


6792 


-16.36 ±0.03 


2.36 ±0.08 


9.86 ±0.32 


0.60 ±0.47 


2.63 ±1. 


.21 


0.98 






o.i z 


8111 


-16.67 ±0.03 


2.98 ±0.08 


9.24 ±0.24 


-0.29 ±0.29 


4.44 ±0 


.79 


1.01 



Note. — / np is the fraction of the luminosity density contributed by the nonparametric fit; in principle / np can be greater than 
unity. The uncertainties are correlated; see Tables 4-8 for the correlation matrix of the uncertainties. See Table 9 for the correlation 
matrix of the luminosity densities and the luminosity evolution parameters of all the bands. 
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Tablc 4. 01 u band uncertainty correlation matrix 





SjM 


SQ 


6P 


(50* 


(5M* 


5a 


SjM 


1.000 


0.949 


-0.934 


-0.938 


-0.091 


0.746 


SQ 


0.949 


1.000 


-0.955 


-0.834 


0.042 


0.823 


SP 


-0.934 


-0.955 


1.000 


0.837 


-0.013 


-0.802 


<50* 


-0.938 


-0.834 


0.837 


1.000 


0.425 


-0.484 


<5M* 


-0.091 


0.042 


-0.013 


0.425 


1.000 


0.560 


6a 




0.823 


-0.802 


-0.484 


0.560 


1.000 



Note. — Correlation matrix calculated from 30 jackknife resam- 
plings of the data, as described in the text. This correlation matrix 
determined from the fl — 0.3, Qa = 0.7 cosmology. 



Table 5. g band uncertainty correlation matrix 





5jM 


SQ 


SP 


(50* 


<5M* 


6a 




1.000 


0.930 


-0.885 


-0.701 


-0.093 


0.493 


6Q 


0.930 


1.000 


-0.949 


-0.494 


0.130 


0.584 


5P 


-0.885 


-0.949 


1.000 


0.447 


-0.144 


-0.656 


<50* 


-0.701 


-0.494 


0.447 


1.000 


0.766 


0.219 


(5M* 


-0.093 


0.130 


HI II 


0.766 


1.000 


0.760 


5a 


0.493 


0.584 


-0.656 


0.219 


0.760 


1.000 



Note. — Correlation matrix calculated from 30 jackknife resam- 
plings of the data, as described in the text. This correlation matrix 
determined from the f2o = 0.3, Oa = 0.7 cosmology. 
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Tabic 6. V band uncertainty correlation matrix 





SjM 


SQ 


5P 




6M. 


5a 


SjM 


1.000 


0.724 


-0.295 


-0.644 


0.325 


0.053 


SQ 


0.724 


1.000 


-0.849 


-0.285 


0.460 


0.222 


SP 


-0.295 


-0.849 


1.000 


-0.048 


-0.364 


-0.296 


<50* 


-0.644 


-0.285 


-0.048 


1.000 


0.498 


0.654 


5M» 


0.325 


0.460 


-0.364 


0.498 


1.000 


0.866 


5a 


0.053 


0.222 


-0.296 


0.654 


0.866 


1.000 



Note. — Correlation matrix calculated from 30 jackknife resam- 
plings of the data, as described in the text. This correlation matrix 
determined from the fl — 0.3, Qa = 0.7 cosmology. 



Tabic 7. H band uncertainty correlation matrix 





5jM 


SQ 


SP 


54>* 


6M* 


5a 




1.000 


0.889 


-0.763 


-0.195 


0.495 


0.568 


6Q 


0.889 


1.000 


-0.950 


0.021 


0.599 


0.607 


5P 


-0.763 


-0.950 


1.000 


-0.131 


-0.574 


-0.654 


f>0* 


-0.195 


0.021 


-0.131 


1.000 


0.735 


0.616 


5M* 


0.495 


0.599 


-0.574 


0.735 


1.000 


0.905 


5a 


0.568 


0.607 


-0.654 


0.616 


0.905 


1.000 



Note. — Correlation matrix calculated from 30 jackknife resam- 
plings of the data, as described in the text. This correlation matrix 
determined from the f2o = 0.3, Oa = 0.7 cosmology. 
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Table 8. z band uncertainty correlation matrix 





5jM 


5Q 


5P 


5</>* 


5AU 


5a 




1.000 


0.328 




-0.312 


0.086 


0.249 


SQ 


0.328 


1.000 


-0.908 


-0.265 


0.085 


0.028 


5P 


-0.110 


-0.908 


1.000 


0.171 


-0.008 


-0.053 


5(f)* 


-0.312 


-0.265 


0.171 


1.000 


0.860 


0.760 


5M* 


0.086 


0.085 


-0.008 


0.860 


1.000 


0.885 


5a 


0.249 


0.028 


-0.053 


0.760 


0.885 


1.000 



Note. — Correlation matrix calculated from 30 jackknife resam- 
plings of the data, as described in the text. This correlation matrix 
determined from the fio = 0.3, Oa = 0.7 cosmology. 



Tabic 9. Uncertainty correlation matrix between luminosity density and evolution fits to all SDSS bands 







a 




SQo.iu 


Sjo.i g 


SQO -lg 




SQa.iy. 


Sjo-H 


SQo.ii 




SQo.i z 


5 jo. 


hi 


0.14364 


1.00000 


0.94949 


0.57194 


0.54212 


0.44706 


0.48121 


0.33898 


0.34904 




0.40158 


SQo. 


hi 


0.88265 


0.94949 


1.00000 


0.58686 


0.59719 


0.43549 


0.41320 


0.26201 


0.28808 


0.10207 


0.35697 


5 jo 


x 9 


0.03854 


0.57194 


0.58686 


1.00000 


0.92950 


0.57137 


0.51576 


0.63139 


0.60449 


0.43418 


0.71319 


SQo 


x 9 


0.51236 


0.54212 


0.59719 


0.92950 


1.00000 


0.63424 


0.54834 


0.62528 


0.64712 


0.30430 


0.71878 


djo 


V 


0.02843 


0.44706 


0.43549 


0.57137 


0.63424 


1.00000 


0.72439 


0.83202 


0.67494 


0.56304 


0.68807 


SQo 


V 


0.29843 


0.48121 


0.41320 


0.51576 


0.54834 


0.72439 


1.00000 


0.80281 


0.86152 


0.23505 


0.82870 


Sjo 


■H 


0.03295 


0.33898 


0.26201 


0.63139 


0.62528 


0.83202 


0.80281 


1.00000 


0.88864 


0.59246 


0.86096 


SQo 


■H 


0.43340 


0.34904 


0.28808 


0.60449 


0.64712 


0.67494 


0.86152 


0.88864 


1.00000 




0.93603 


Sjo. 


l z 


0.02110 




0.10207 


0.43418 


0.30430 


0.56304 


0.23505 


0.59246 


0.21461 


1.00000 


0.32770 


SQo. 


!* 


0.29156 


0.40158 


0.35697 




0.71878 


0.68807 


0.82870 


0.86096 


0.93603 


0.32770 


1.00000 



Note. — Correlation matrix calculated from 30 jackknife resamplings of the data, as described in the text. This correlation matrix 
determined from the (£!o = 0.3, Ha = 0.7) cosmology. We keep enough significant figures that no clement of the inverse of this 
correlation matrix differs by more than a couple of percent from that calculated using the machine precision result, purely to avoid 
confusion between results wc might obtain from a machine precision version of this matrix and those others might obtain from this 
table. 



Table 10. Luminosity density if-corrected to various bandpasses 



j + 2.51og 10 /i 



Band 


Aoff (A) 


Atuab 


This paper 


This paper SDSS Comm. Data 


2dFGRS 


2dFGRS 








z = 


z = 0.1 


z = 


z = 


z = 0.1 


o.o u 


3538 


0.00 


-13.99 


-14.35 


-15.21 






o.o g 


4664 


0.00 


-15.27 


-15.46 


-16.05 






o.o r 


6154 


0.00 


-15.90 


-16.06 


-16.41 






0.0, 


7471 


0.00 


-16.29 


-16.41 


-16.74 






o.o z 


8922 


0.00 


-16.74 


-16.79 


-17.02 






o.qg 


3974 


-0.12 


-14.98 


-15.17 








O.Oy 


4980 


0.02 


-15.69 


-15.85 








a.o R 


5906 


0.21 


-16.23 


-16.39 








O.Oj 


7291 


0.45 


-16.90 


-16.99 








°\ 


4141 


-0.08 


-15.18 


-15.36 


-15.97 


-15.35 


-15.45 



Note. — We have taken the luminosity densities in Table 3 for the (Qo = 0.3, Oa) cosmology 
and applied the methods of Blanton et al. (2002a) in order to evaluate the luminosity density in a 
number of other bandpasses. We also show 2dFGRS results from Norberg et al. (2002) and old SDSS 
results from Blanton et al. (2001). We have inferred the z = 0.1 value of the luminosity density from 
Norberg et al. (2002) based on the evolutionary corrections in their Figure 8. BVRI magnitudes 
assume the Bessell (1990) response curves. We set bj = B — 0.28(£> — V) by definition (except in the 
value from Blanton et al. 2001, where we copy the number directly from their table). ArriAB is the 
offset to apply to translate the listed magnitude into an AB magnitude (using the Hayes 1985 Vega 
spectrum). 
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Large-Scale Structure sample 10 

Fig. 1. — Equatorial distribution of right ascension and redshift for Main Sample galaxies within 6° of the 
Equator in sample 10 of the SDSS. The half of the survey in the Galactic South (on the right) appears less 
dense than the Galactic North (on the left) simply because the imaging near the equator extends less in the 
declination direction in the south than in the north. 
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Fig. 2. — Distribution in Galactic coordinates of spectroscopically confirmed galaxies in the SDSS Large- 
Scale Structure sample 10. The effective area covered (the integral of the SDSS sampling fraction over the 
covered area) is about 1844 deg 2 . 
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Wavelength A (A) 



Fig. 3. — Demonstration of the differences between the unshifted SDSS filter system ( Q u, g, ( V, i, z) 
in the top panel and the SDSS filter system shifted by 0.1 (°'V 01 g, ° V, 0m H, 01 z) in the bottom panel. 
Shown for comparison is a 4 Gyr-old instantaneous burst population from an update of the Bruzual A. & 
Chariot (1993) stellar population synthesis models. The if -corrections between the magnitudes of a galaxy in 
the unshifted SDSS system observed at redshift z = 0.1 and the magnitudes of that galaxy in the 0.1-shifted 
SDSS system observed at redshift z — are independent of the galaxy's spectral energy distribution (and 
for AB magnitudes are equal to — 2.51og 10 (l + 0.1) for all bands; Blanton et al. 2002a). This independence 
on spectral type makes the 0.1-shifted system a more appropriate system in which to express SDSS results, 
for which the median redshift is near redshift z = 0.1. 
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0.05 0.10 0.15 0.20 

Redshift z 

Fig. 4. — if-corrections for the twelve SED types used for this paper, for each band. Because we are 
i^T-correcting to bands shifted to z = 0.1, all galaxies have the same if-correction (— 2.51og 10 (l + 0.1)) at 
z = 0.1. For this reason, choosing this set of bandpasses minimizes our uncertainties in the luminosity 
density at z = 0.1. 
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-24 -22 -20 -18 

M 0Ar - 5 log 10 /i 

Fig. 5. — Luminosity function in the ° V band. The thick solid line is the luminosity function fit; the thin solid 
lines are the individual gaussians of which it is composed. The grey region around the luminosity function 
fit represents the lcr uncertainties around the line; naturally, these uncertainties are highly correlated with 
each other. The dashed line is the Schechter function fit to the result. The luminosity density the evolution 
parameters, and the parameters of the Schechter function are listed in the figure. 




Fig. 6. — Same as Figure 5, for the 01 u, 01 g, °'H, and 01 z bands. 
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Thick Histogram: Data (147986 galaxies) 
Thin Histogram: Model (150316 galaxies). 
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Fig. 7. — Redshift distribution of the ° V band sample, for each quartile (weighted by number) in absolute 
magnitude. The thick line represents the data; the thin line is a Monte Carlo representation of the model, 
including the selection effects in the survey. In this figure and in the following Figures 8-12, the model 
is a decent representation of the data, but not a perfect one; much of the difference is likely to be due to 
large-scale structure, but it is possible that further complications of our evolution model or our error model 
might be necessary to fully reproduce the data. 
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Fig. 8. — Absolute magnitude distrirjuifeAxfai the °Q banc?, Sot Several slices of rcdshift. The thick line 
reperesents the data; the thin line is a Monte Carlo representation of the model, including the selection 
effects in the survey. 
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Fig. 9.— Same as Figure 8, for the ° ^>M(g ~ ^ ^ ° § 1 0^ 
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Fig. 10. — Same as Figure 8, for the 
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Fig. 11.— Same as Figure 8, for the °4^tPaM. 
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Fig. 12.— Same as Figure 8, for the ( M%Aa&. ~ ^ ^°§>10^ 1 
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Fig. 13. — Differences Am = ms -mp between Sersic and Petrosian magnitudes for each band for a volume- 
limited sample with —23 < Mo.^. < —21 and 0.02 < z < 0.15. The left column shows the differences as a 
function of Sersic index n. As expected, for galaxies at high Sersic index (close to the de Vaucouleurs value 
n = 4) the Petrosian magnitudes are an underestimate relative to Sersic magnitudes. For r, i, and z, the 
differences between the two remain small at low Sersic index (near the exponential value n = 1). The right 
column shows the differences as a function of redshift z. A linear regression is shown as a grey line, along 
with the parameters associated with the best-fit regression. The slopes are generally insignificant compared 
to our uncertainties in the evolution parameter Q. 
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Luminosity Density vs Wavelength 

O SDSS at 2=0.1 
A 2dF at 2=0.1 
X 2MASS at z=0.1 
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Fig. 14. — Luminosity density and its evolution as a function of wavelength, for the SDSS (diamonds; this 
paper), the 2dFGRS (triangle; Norberg et al. 2002), and 2MASS (crosses, Cole et al. 2001). The 2dFGRS 
has been evolution-corrected to z = 0.1 using Qi, j = 1 (the effective value used by Norberg et al. 2002) and 
2MASS has been evolution-corrected to z = 0.1 using Qj = Qk = 1- 



